Number of branches in diffusion-limited aggregates: The skeleton 
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Abstract 



We develop the skeleton algorithm to define the number of main branches 
A*";, of diffusion-limited aggregation (DLA) clusters. The skeleton algorithm 



^ ' provides a systematic way to remove dangling side branches of the DLA clus- 



ter and has successfully been applied to study the ramification properties 
of percolation. We study the skeleton of comparatively large (~ 10^ sites) 
off-lattice DLA clusters in two, three and four spatial dimensions. We find 
that initially with increasing distance from the cluster seed the number of 
branches increases in all dimensions. In two dimensions, the increase in the 
number of branches levels off at larger distances, indicating a fixed number 
of = 7.5 lb 1.5 main branches of DLA. In contrast, in three and four di- 
mensions, the skeleton continues to ramify strongly as one proceeds from the 
cluster center outward, and we find no indication of a constant number of 



main branches. Likewise, we find no indication for a fixed Ni, in a study of 

DLA on the Cayley tree. In two dimensions, we find strong corrections to 

scaUng of logarithmic character, which can help to explain recently reported 

deviations from self-similar behavior. 
61.50.Cj, 05.40.+j, 64.60.Ak, 81.10.Jt 
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I. INTRODUCTION 



The fractal nature of diffusion-limited aggregation has been established by a number of 
studies, and a theoretical framework is being developed In the initial paper written by 
Witten and Sander [jl| , power-law behavior of the two-point correlation functions was found. 
Later, Meakin ||^ observed power-law scaling of the radius of gyration of DLA in two (2D), 
three (3D), and four dimensions (4D). Enhanced on- and off-lattice algorithms facilitated 
studies of larger aggregates (Meakin [Q, Ball and Brady |^) and a striking influence of 
the lattice structure on the structure of the aggregates was found Advanced off-lattice 
simulations by Ossadnik grew clusters up to M = 5 x 10^. All these studies indicate or 
confirm a fractal dimension of df = 1.715 ± 0.004 as reported by Tolman and Meakin [^f. 
Apart from the scaling of the radius of gyration with cluster mass, fractality is revealed in 
studies of the branching structure of the aggregates by, e.g., Alstr0m et al. 0, Hinrichsen 
et al. |10|, Ossadnik [|ll[|, and more recently Yekutieli et al. |[12|. Self-similarity has also 
been found in a study by Argoul et al. calculating the moments of the mass distribution 
of on-lattice clusters [0, by Barabasi and Vicsek's mapping of the cluster perimeter to a 
self-affine surface ||14[, in investigations of the fjord geometry of DLA [|l^,|l^, and many 



more (for reviews, see, e.g, p]7|-p2|). 

However, almost as old as the model itself are reports about deviations from a simple 
statistically self-similar structure. In Refs. [^p3| lattice anisotropy and deposition habit 
effects on the structure of lattice clusters were studied, Meakin and Vicsek [^] found that 
the tangential and radial two-point-correlation functions in DLA display different exponents. 



In recent works however, Hegger and Grassberger ||25| point out that off-lattice DLA in 
strip geometry displays a slow approach to local isotropy. Differences arise in the fractal 



dimensions measured in circular or strip geometry ||26|. Simple self-similar scaling is not 
sufficient to fully characterize the mass distribution within a DLA cluster, as recognized by 



Meakin and Havlin Vicsek, Family and Meakin ||2^, Amitrano, Coniglio, Meakin, and 
Zanetti and Mandelbrot [^. Mandelbrot examined the lacunarity pT|j3^ of off-lattice 
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DLA and found a trend towards increasing compactness for increasing systems size which is 
not reflected in the measurements of the radius of gyration of the clusters. 



In our study we will apply the "skeleton" algorithm [^,0] to calculate the number of 
main branches Ni, in off-lattice DLA (Sec. |I|). We obtain the skeleton of a DLA cluster by 
removal of all the small dangling side branches of the cluster (we define the skeleton in the 
following section). Only the main branches that reach out to a certain fraction of the cluster 
radius remain. For a self-similar, scale-invariant structure we expect that the large-scale 
properties, like the number of main branches of the cluster, are not altered under rescaling. 
Conversely, we take an increase or decrease of the number of main branches with increasing 
cluster size to suggest that the structure is not self-similar. 

In 2D, we find that significant finite size effects of logarithmic nature are present which 
may be related to some of the observations reporting deviations from simple self-similar 
behavior. Asymptotically however, the number of main branches of 2D DLA approaches 
a constant number (Sec. |III| ). In contrast to the behavior in 2D, the number of branches 
of 3D and 4D DLA increases with cluster mass, indicating that these clusters may not be 
self-similar (Sec. fV\ ). 

II. THE SKELETON 

DLA has a loopless tree structure [^. This observation is apparent in simulations of 
off-lattice aggregates ||36| , |37[| where an incoming particle is added to the cluster when its 



distance to the cluster is below a specific sticking distance. Its "parent" particle is the one 
that it was closest to at the moment of incorporation into the cluster. The child-parent 
relationship allows us to uniquely assign a generation number or "chemical distance from 
the seed" to every cluster particle. To this end, we assign to the seed particle of the DLA 
cluster the number i = 0. The children of the seed are assigned the chemical distance 
i = 1] children of £ = 1 particles have chemical distance i = 2. In general, each particle 
inherits the chemical distance of its parent and increases this number by one to find its own. 
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The chemical distance turns out to be a very useful quantity in studies of the branching 
properties of DLA |P-p!l|] and we will use it here to define the skeleton of DLA — closely 



following Ref. |3|. 

Roughly speaking, we obtain the skeleton by removal of all the small dangling side 
branches of the cluster. To this end we first identify all the sites that are tips of branches 
— those sites to which no other particle has been attached during the growth process. The 
^ value of each tip £tip is then passed back to its parent and grandparent and so forth until 
we hit an ancestor site which has more than one child. In this case, only the largest of the 
£tip values of the children is retained. We iterate this process until all cluster sites "know" 
their £tip value. 

Now let us choose an arbitrary value ic less than the largest chemical distance A of any 
site from the seed. Given "we define as the skeleton those sites for which both £tip > and 
i <lc holds. Thus, side branches which have not grown out to at least Ic do not contribute 
to the skeleton. 

In Fig. |T], we display the skeletons of two randomly selected growing DLA clusters at 
mass M = 5 000, 50 000, 500 000. In this figure, 4 is chosen to be half the "chemical radius" 
A/2 of the cluster. Note that the termination points of the skeleton are located almost on 
a circle, although we use the chemical and not the Euclidean distance from the seed to 
define the skeleton. This indicates that DLA grows radially outward without forming loops. 
Indeed, it has been found that the fractal dimension rfmin of the path connecting the cluster 
seed to a given site is equal to 1 Here, rfmin is defined by £ ~ r'^™'". In other words, for 



cluster sites at a given distance r from the cluster seed, the average value (£) is proportional 
to r, (£) ~ r. For the remainder of this study, we will consider i and r as equivalent. 

One physical interpretation of the skeleton in DLA can be obtained if we consider the 
aggregate as a conductor situated between the grounded seed and a circular electrode (a 
sphere or hypersphere in 3D or 4D, respectively) of "radius" ic- Then the skeleton is the 
collection of paths that contribute to the current through the aggregate. 

The skeleton as defined above has some desirable and some less desirable features. First, 



its definition is certainly simple, which is desirable. Second, we see in Fig. |T] that in an 
intermediate range of ^ the number of branches in the skeleton is small, which is another 
desirable feature since we claim that the skeleton identifies the main branches of the DLA 
cluster. However, we see that close to £c a lot of branches appear. Their presence ultimately 
reflects our ignorance as to whether a specific side branch will keep growing or die out. Thus, 
these branches are a very "physical" property, but, nevertheless, an undesirable feature. We 
will address this problem by keeping £c a variable and systematically test the dependence of 
the skeleton calculations on this parameter. 

We define A = £c/A and consider ensembles of DLA clusters for which the the value A 
is fixed. For such an ensemble, we define -B(A, A; t) as the average number of sites in the 
skeleton at chemical distance d from the cluster seed. Note that i?(A,A;£) is defined only 
for £ < 4. 

Some general conclusions about properties of the DLA skeleton can be drawn from the 
definition and are visible in Fig. |1|. 

(i) For £ = 0, only the cluster seed contributes, thus B{X,A;i = 0) = 1 independent of 
spatial dimension d and cluster size A. 

(ii) Since i ^ r, and since DLA is essentially loop less so that overhangs can be neglected, 
we can determine the value B{X,A;i = ic) from the number of intersection points of 
a "cookie cutter" of radius corresponding to ic with the DLA cluster. As the analog 
of a cookie cutter in 3D we consider a sphere, and in 4D a hypersphere. The result 
is B{\,A;£c) ~ ^c'~^ where df denotes the fractal dimension of DLA, which depends 
on the embedding dimension d. The exponent df — 1 follows from the fact that the 
codimension of the intersection of two sets is obtained as the sum of the codimensions 
of the two intersecting sets. 

(iii) The number of sites in the skeleton increases monotonically with i, i.e., -B(A, A; i') > 
B{X,A;i) if i' > i. Qualitatively we observe from Fig. |I| that the increase is slow in 
the vicinity of the cluster seed, but is followed by a precipitous increase close to ic- 
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(iv) For fixed £ and cluster "radius" A, the number of sites in tlie skeleton decreases mono- 
tonically when A is increased, because some branches that contribute to the skeleton 
at some value of ^c fail to reach larger values. In the extreme case, when A = 1, we 
obtain 5(1, A; €) = 1, independent of A and i. 

In the study of the skeleton of percolation clusters a slow increase of the number of 
arms for small i was found to be due to finite size corrections. Taking this into account, one 
obtains for small i a constant number of branches. In analogy to the percolation example, 
we use finite-size scaling techniques to extract information about DLA from our 2D, 3D and 
4D skeleton calculations, which we present in the following sections. The DLA results in 
finite dimensions will be compared to the skeleton of a random Cayley tree model, which 
can be considered to be the limit of infinite dimensionality of DLA [p^j . 



III. SKELETON OF 2D DLA 

In Fig. I we display B{\ = 0.5, A; £) as a function of £/4 for A = 100, 316, 1 000, 3 162, 
and 5 000. For A < 1000, we obtain our data from 500 off-lattice clusters, for A > 3162, 
114 clusters are averaged. The clusters are grown stepwise and growth is interrupted for the 
analysis when A equals one of the above listed values. For A = 5 000, the cluster mass M is 
typically M = 1 - 1.5 x 10^ sites. 

As a function of i/ic we find that after an initial short transient the central part of 
the data displays approximately linear behavior. The linear region widens when the cluster 
size A increases. The approximate linear behavior terminates when i/ic approaches ~ 0.3. 
Beyond i/ic ~ 0.3, B{0.5, A; i) increases sharply towards the final value ~ it^ ^. However, 
the linear region of the data tends to become flatter as we proceed to larger A. Since -B(A, A) 
is monotonically increasing, the slope a;(A, A) in the linear region cannot be negative. Thus, 
asymptotically q;(A, A) can either approach or assume a finite value > 0. We will refer to 
the restricted range of approximately linear behavior as the "scaling regime" of the skeleton. 

In order to distinguish between these two possibilities, we examine the dependence of 
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a{X, A) on system size, here characterized by A. To this end, we determine the slope a(A, A) 
of a hnear least squares fit. 

In B{X, A; i) ^ a{X, A) ln(|) + (3{X, A), (1) 

to ln_B(A,A;£) vs \n{i/ic) for data points satisfying i > 10 and i/£c < 0.3. We use the as- 
sumption, familiar from the estimate of critical exponents from finite system size calculations 
(see, e.g. |3^), that a{X,A) may be written as 

a(A,A) = aoo(A)+/(A,A) + .... (2) 

Here, the leading correction term /(A, A) to the asymptotic value aoo(A) is typically either 
(i) a power law ~ A^'^ or (ii) a logarithmic term, e.g., ~ 1/lnA. Both vanish for A — > oo. 
The dots indicate corrections that decrease faster than /(A, A) as A — >■ oo. 

In Fig. 1^ we plot a;(A;A) vs 1/lnA for different values of A. From the almost linear 
behavior of a;(A;A) for large A we conclude that the leading order correction term can be 
approximated by the reciprocal of In A (ii), i.e., 

/(A,A)=A(A)/lnA. (3) 

Here, A{X) denotes an amplitude which depends on A as evident from Fig. |^a. However, our 
calculations do not allow us to exclude the possibility that /(A, A) vanishes in the power-law 
fashion (i), although, in this case the characteristic power 7 must be small, 7 < 0.2. Larger 
values of 7 would lead to significant curvature in the large A behavior of a;(A,A) [cf. Fig. 
|b]. 

From figure |^ we conclude that aoo(A) = for 0.3 < A < 0.8 by extrapolating the 
trend displayed by a(A,A) for the A range accessed by our simulation. Since for A — 1 , 
B{X, A;i) = l for all i, we expect this result to also hold if A > 0.8. For small A, the skeleton 



probes the frozen region of the DLA cluster which effectively stopped growing so that 
our results should be less effected by possible slow structural changes in the growth zone of 
the cluster |]TB|,^^. Therefore, although the fit region of B{X,A;i) becomes narrower, we 
assume that aoo(A) = extrapolates into the small A regime. 
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The implications of aoo(A) = are intriguing. Asymptotically, we can assign a fixed 
number of main branches A'^f, to 2D DLA growth. We deduce the actual number of these 
branches from the behavior of the intercepts l3{X; A) of the fit ([1|) in the "fiat" region of the 
skeleton. In Fig. |^c, /5(A;A) is plotted as a function of 1/lnA for several A. The value of 
(3{X; A) slightly increases as the cluster size grows — corresponding to more main branches 
in larger systems. However, the increase levels off at about /?oo(A) = 2 ± 0.2, so that on 
average A";, = exp[/?oo(A)] = 7.5 ± 1.5 main branches result. 



IV. DLA IN HIGHER DIMENSIONS: 3 AND 4 

We also perform simulations of off-lattice DLA in 3D and 4D. Our findings are quite 
different from the 2D case, a result not surprising in the light of earlier work suggesting that 



2D is a special case [^. In particular, our calculations suggest that in contrast to the 2D 
case, where aoo(A) = 0, in the 3D and 4D case aoo(A) > 0. We discuss in the following the 
numerical data supporting this result. 



1. Skeleton of 3D DLA 

In Fig. ^ we display B{X = 0.5, A; £) for 3D DLA. Different curves are characterized by 
different values of A. The largest value of A is 560, which corresponds to DLA clusters with 
M ^ 10^ particles. The skeletons of 450 (A < 200) and 75 (A > 200) clusters have been 
averaged to obtain the displayed data. 

The tendency for the central part of the plot to become fiatter as A increases is much 
less pronounced than in 2D. This observation is supported by Fig. which shows the 
A dependence of q;(A;A). We obtain the slope q;(A;A) from a linear least square fit of 
In B{X, A] i) in the region characterized by ^ > 10 and i/i^ < 0.3. Like in 2D, also in 3D 
the correction term is ~ 1/lnA (^, and likewise we cannot exclude that the correction is 
of power-law type with a small exponent. However, unlike in 2D, in the 3D case the slopes 
corresponding to the small A values associated with the "frozen" region of the cluster a (A; A) 
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extrapolate to a non-zero value, aoo(A) ~ 0.2. This value indicates that in 3D the number of 
branches of the skeleton increases monotonically in power-law fashion as a function of 1/ 1^. 
even in the limit of infinite cluster size. 

2. Skeleton of 4D DLA 

The skeletons of 50 off-lattice DLA clusters were averaged to obtain the data displayed 
in Fig. ^. The clusters with A = 190 contain up to 500 000 particles. Similar to the 
3D case, the slopes in the central part of the i?(A,A;£) data do not tend to become flat 
with increasing mass and this behavior is reflected in the A dependence of a(A;A), which 
we display in Figs, ^p, and b. The quality of the 4D calculations is much poorer than that 
of the 3D and 2D calculations. A finite value aoo(A) is consistent with the two finite size 
scaling plots Figs. |^ and b. As in 3D, such a finite value indicates a power-law branching 
of the skeleton in the scaling regime as the cluster mass goes to infinity. 

Moreover, if we consider the intersections /5(A,A) of the straight line fits to the scaling 
region we observe that they may diverge as A oo. Such behavior would imply that 
i?(A, A; £) may not converge to a limit as A — oo. Thus we expect rather strong deviations 
from self-similar behavior in 4D and for A ^ oo no self-similar limit as in the 2D case may 
exist. 

V. CAYLEY TREE MODEL 

The realization of the DLA growth model on a Cayley tree is considered to be the mean- 
field or infinite dimensionality limit of DLA. DLA on the Cayley tree has been formulated 
as follows |3^. To the cluster seed (shell 0), we links z surface sites (shell 1) — z > 2 is 
the functionality of the model. In each growth step, we select (occupy) one of the surface 
sites randomly. We then create z — 1 new (empty) surface sites and link them to the just 
occupied site. The shell number of the new empty sites is equal to the shell number of the 
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occupied parent incremented by 1. Note that all occupied sites have z links to neighbor sites 
and that all empty sites have exactly 1 neighbor site. 

Since new cluster sites are randomly picked from the set of all surface sites of the growing 
tree, we effectively disregard screening of the interior parts of the tree, which in finite 
dimensions is the essential mechanism to produce low dimensional fractal structures [03 



We have numerically grown clusters of this tree model with functionality z = 3, 4 and 5 as 
large as 32 shells. As a function of the shell number, the number of occupied sites grows 
exponentially ^Sf. Similar to the case of 3D and 4D DLA, the skeleton shows no sign of 



saturation as a function of i/ic- The increase seems to have even exponential character, 
with ln_B(A, A = 32; i) oc i/ic, and prefactors that depend on z and A. 

This behavior is consistent with the situation in 3D and 4D DLA, where the skeleton as 
a function of i/ic also increases, albeit in a power-law fashion. 

VI. DISCUSSION 
A. Previous results on the number of stable branches of DLA 



The number of branches A^;, of DLA has been used as an important morphological char- 
acterization of the cluster In 2D using conformal mapping arguments. Ball \W 



has proposed that the maximum number of stable branches A^max is related to the mass 
scaling dimension D of the cluster by 

(Ar^ax/2-l)(D-l) = L (4) 



For D ^ 1.7 this relation predicts A^max ~ 4.9. Similarly, Turkevich and Scher have 
modeled the the tips of the cluster as the corners of a convex polygonal cluster envelope. 
The number A^ of branches then determines the tip angles and those in turn the mass scaling 
dimension, namely, D — 1 = N/{N + 2). If one replaces D — 1 from Eq. and assumes 
A^ = A'max, we obtain A'max ~ 4.8, in agreement with the previous prediction. 
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Both numbers seem to be significantly smaller than our finding Nb7.5 ±1.5 in 2D.Q 
It is not clear why Nf, should be different from A^max following from Eq. On the one 
hand, if the actual number of branches is larger than the maximum number of stable ones 
then the competition among branches for the incoming flux of random walkers will cause 
branches to die (cf. If, on the other hand, the number of branches is too small, then 

new branches will be created because the fjords of the aggregate are not sufficiently screened 
to suppress growth. If, in fact, A^max and Nb do not coincide, the reason could be that large 
branches may finally split, and one of the now too numerous arms will finally die out, but 
slowly enough to increase the number of branches apparent in the skeleton. 



B. Relation to "lacunarity" measurements on DLA 



Recently, Mandelbrot has investigated the lacunarity of off-lattice DLA |^. To calculate 
the lacunarity of a cluster of span A, we (i) select the frozen interior region of the cluster (of 
span kA, k < 1) and rescale all the particles coordinates into the unit circle, then (ii) center 
solid disks of radius e ^ 1 on the rescaled particle positions and (iii) determine the area 
fraction / of the unit circle covered by the e disks. For a "normal" fractal, one expects / 
to saturate when A increases and in this case to provide a measure for the lacunarity of the 
object. For DLA however, / increases continuously with A. Mandelbrot argues accordingly 
that either DLA displays a "massive transient" or a "limitless drift" towards increasing 
compactness. 

We argue that the logarithmic finite size corrections that we see in 2D in the approach of 
B{X,A;i) to its asymptotics can explain the "massive transient" picture. As we have seen 
above, B{X,A;i) increases with system size corresponding to structurally different clusters 



Concerning the Turkevich-Scher assumption of a convex polygonal cluster envelope one may 
argue that it is possible to conceive cluster envelopes with more acute angles and a larger number 
of arms — possibly around 7. 
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with more main branches Ni,. Consequently, the area fraction / increases. 

Thus, off-lattice DLA shows a very slow structural change towards its asymptotic struc- 
ture, similar to its on-lattice siblings. 



C. Why is 2D DLA different from higher dimensional DLA? 

The remaining open question is why on the one hand aoo(A) — > in 2D, but, on the other 
hand, aoo(A) > in 3D and 4D and in the latter case possibly even a divergence of /5(A, A) 
with A. Let us first note, that 2D is known to be a "critical dimension" for radial DLA with 



respect to the scaling properties of the growth probabilities of the surface sites P2| , |51| , |52 
We find two more indications supporting our result, one from a simple consideration of 
the cluster density in different spatial dimensionalities, the other from a comparison of the 
correlated DLA tree to a random tree model, for which the scaling behavior of the skeleton 



is known p4|. 



1. Cluster "density" and screening 

We have already noted above that the number of branches in the skeleton at its termi- 
nation points at is ~ £c'' ^. Thus, for given £ = there are many more branches present 
in the 3D and 4D skeleton than in the 2D skeleton. In order that these branches contribute 
to the skeleton in the scaling region they must originate deep in the cluster interior. 

With increasing dimension, the screening of the cluster interior from incoming particles 
decreases. This decrease manifests itself for example in the strong semi-exponential screening 
of the cluster fjords in 2D |[53|-|57| in contrast to the power-law type screening in 3D or 



in the screening free Cayley-tree model — where we also observe aoo(A) > 0. Qualitatively, 
the reduction of screening can also be seen in the mean-field expression for the particle 
penetration depth ~ p^^ ~ A'^^'^f . The penetration depth increases with d, because df 
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approaches d — 1 from above.0 

Therefore we conclude that the weaker screening in dimensions > 2 is responsible for 
the further increase in the number of branches of the skeleton both as a function of i/ic and 
as a function of cluster size A. 

2. Comparison to a random tree model 
DLA has the structure of a tree. Thus, we find it instructive to compare our findings 



for the off-lattice DLA skeleton to results for the skeleton of a random tree model |^ 
This tree model is defined on a 2D square lattice and characterized by a tunable "intrinsic" 
dimension de, which determines how the mass M{i) of the tree object depends on the 
chemical distance i from the cluster seed, i.e., M{i) oc f*^^. To this end, we occupy randomly 
exactly N{£) = sites in the ith chemical shell (or fewer if necessary) — all the other 
sites are blocked. Only those sites in the ith shell shall be occupied that do not close 
loops, i.e., that are not neighbors to sites in previous shells. Note that apart from possible 
constraints due to the embedding lattice at small i, the topological structure of the clusters 
does not depend on the dimension of the embedding lattice. 

Reference then considers the skeleton 5(1, A; of a tree with intrinsic dimension d^. 



In contrast to DLA, here in general B{l,A;i) > 1, because the last chemical shell contains 
A^(A) = ^"^^"^ ^ 1 sites. In these random tree models there exists a sharp transition value 
dj = 1.65 ±0.05 which separates two regions with different types of behavior. For trees with 
di < d1 the skeleton displays a flat — approximately constant — region as a function of 
like the one observed in the 2D DLA case {non-branching skeleton). However, for larger 
values of di the skeleton displays a well deflned power-law increase as function of £, i.e., 
ttoo(l) > [branching skeleton). 



^ It is however known from simulations that at least in 2D this relation does not hold and the 
penetration depth is ~ A. 
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Thus, if we consider DLA as such a random tree, we expect a branching skeleton for 
spatial dimension d > 2, since = df > ri^. In 2D, df is only slightly larger than d^, or 
even equal to dj within the error bars, which is consistent with the observed non-branching 
skeleton. 

VII. SUMMARY 

We have determined the skeleton of comparatively large (~ 10® sites) off-lattice DLA 
clusters in 2D, 3D and 4D. We find that, asymptotically, in 2D the skeleton of DLA suggests 
a fixed number of Nj, ~ 7.5 ± 1.5 main branches and a self-similar structure. In 3D and 
4D, and possibly all spatial dimensions o? > 4, the DLA skeleton is a ramified object, which 
displays branching over the whole range of £- values for which it is defined. 

For all dimensions we find strong finite-size effects corresponding to a slow change in 
the structure of DLA as the cluster size increases. The presence of strong corrections to 
scaling of logarithmic character in 2D DLA is in agreement with findings of deviations from 
self-similar behavior reported elsewhere. 

Moreover, in 4D, it is possible that the structure of the aggregate keeps changing even 
asymptotically, such that there is no self-similar limit. 
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FIGURES 

FIG. 1. Skeleton for 2 growing DLA clusters: (a-c) first cluster, (d-f) second cluster. Growth 
of the cluster has been interrupted at cluster masses M = 5 000 (a,d), 50 000 (b,e) or 500 000 (c,f). 
The determination of the skeleton is based on a value of ^c/A = 0.5. The skeleton is then rescaled 
so that independent of cluster mass the same size results. 

FIG. 2. Log-log plot of the 2D skeleton B{\ = 0.5, A; £) vs £/£c for different A values as 
indicated in the legend. For A < 1 300 data are averaged over 850 off-lattice DLA, for A > 1 300 
averages over 114 clusters are taken. For comparison, we have also averaged data in the con- 
stant M ensemble, instead of the constant A ensemble used here, and obtain similar results (not 
shown). Note that in an intermediate region of the plot {(/£c < 0.3) a slope can be associated 
with B{\ = 0.5, A; ^) vs. £/£c- As A increases, this slope becomes increasingly smaller, indicating 
that the asymptotic curve may be fiat in this "scaling" region. Inset: The number of branches in 
the skeleton at ^ = £c, B{X,A;£c) is equal to the number of intersection points of a cookie cutter 
of radius £c in £ space with the cluster, i.e., B{X,A;£c) ^ A'^f~^. The solid line is a guide to the 
eye and indicates a slope of — 1 = 0.70. Different symbols refer to different A = 0.8(0), 0-6(n), 
0.5(A), 0.4(v), 0.2(.). 

FIG. 3. Finite size behavior of the slopes a(A, A) of B{X = 0.5, A; ^) vs £/ic calculated from 
data in the range ^ > 10 and £/£c < 0.3. The slopes a(A, A) for different A as indicated in the 
legend of the figure are plotted vs (a) 1/lnA and (b) 1/A. In (a) the slopes lie on asymptotically 
straight lines indicating a finite size correction term ~ 1/lnA. This behavior should be contrasted 
to the calculations in (b) where the correction term is assumed to be proportional to 1/L and 
pronounced curvature is displayed as A increases, (c) Intercepts /3(A, A) of the straight line fit with 
the ordinate vs 1/lnA. If aoo(A) = 0, then /3oo(A) denotes the logarithm of the number of arms 
Nh in the "flat" region of the skeleton. From jS^ « 2 ± 0.2 we conclude that = 7.5 ± 1.5. 
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FIG. 4. Log-log plot of the 3D skeleton B{X = 0.5, A; £) vs £/4 for different A values as 
indicated in the legend. We take averages over 75 (A > 200) and 450 (A < 200) clusters. Note that 
the tendency of the power-law region in the plot to become flatter with increasing A is much less 
pronounced than in the 2D case |^ and is in fact rather indicating a non-zero aoo(A) > 0. Inset: 
Number of branches in the skeleton at £ = £c- We expect B{X,A;£c) ~ A'^^'^"'^. The solid line is a 
guide to the eye and indicates a slope of 1.5 ~ d/ — 1 in 3D. Different symbols refer to different A 
and are the same as in Fig. 

FIG. 5. Finite size behavior of the slopes a(A, A) of B{X = 0.5, A; ^) vs i/ic calculated from 
data in the range i > 10 and i/ic < 0.3 for 3D DLA. The slopes a(A, A) for different A with 
symbols as in Fig. ^ are plotted vs (a) 1/lnA and (b) 1/A. (c) Intercepts /3(A, A) of the straight 
line fit with the ordinate vs 1/ In A. 

FIG. 6. Log-log plot of the 4D skeleton B{X = 0.5, A; £) vs £/ic for different A values as 
indicated in the legend. As in the 3D case, the central region of the plot does not display a 
tendency to become flat as A increases. However, the system sizes seem still too small to draw a 
final conclusion, in particular if one considers the size dependence of the slopes in the power-law 
region and the comparatively poorer quality of the calculations (see Fig. |^. Inset: Number of 
branches in the skeleton at i = £c- We expect B{X,A;£c) ~ A"'^'^-'^. The solid line is a guide to the 
eye and indicates a slope of 2.2 ~ df — 1 in 4D. Different symbols refer to different A and are the 
same as in Fig. ^. 

FIG. 7. Finite size behavior of the slopes a(A; A) of B{X = 0.5, A; £) vs £/ic calculated from 
data in the range i > 10 and i/ic < 0.3 for 4D DLA. The slopes a{X, A) for different A with 
symbols as in Fig. ^ are plotted vs (a) 1/lnA and (b) 1/A. (c) Intercepts /5(A, A) of the straight 
line fit with the ordinate vs 1/ In A. 
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